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ABSTRACT 

We report new calculations of interstellar 15 N fractionation. Previously, we have shown 
that large enhancements of 15 N/ 14 N can occur in cold, dense gas where CO is frozen out, but 
that the existence of an NH + N channel in the dissociative recombination of N2H + severely 
curtails the fractionation. In the light of recent experimental evidence that this channel is in 
fact negligible, we have reassessed the 15 N chemistry in dense cloud cores. We consider the 
effects of temperatures below 10 K, and of the presence of large amounts of atomic nitrogen. 
We also show how the temporal evolution of gas-phase isotope ratios is preserved as spatial 
heterogeneity in ammonia ice mantles, as monolayers deposited at different times have differ- 
ent isotopic compositions. We demonstrate that the upper layers of this ice may have 15 N/ 14 N 
ratios an order of magnitude larger than the underlying elemental value. Converting our ra- 
tios to (5-values, we obtain S 15 N > 3, 000 %c in the uppermost layer, with values as high as 
10,000 %o in some models. We suggest that this material is the precursor to the 15 N 'hotspots' 
recently discovered in meteorites and IDPs. 
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1 INTRODUCTION 

Laboratory analyses of primitive solar system materials, such as 
meteorites, interplanetary dust particles (IDPs), and cometary dust 
particles returned by the Stardust mission, show anomalous frac- 
tionation in the heavy isotopes of numerous elements relative to 
that expected from the cosmic or solar system values (Clayton & 
Nittler 2004; Lodders & Amari 2005; Alexander et al. 2007; Mc- 
Keegan et al. 2006). In the case of hydrogen and nitrogen, the large 
D/H and 15 N/ 14 N ratios observed in some phases have been at- 
tributed to the survival of D- and J N-enriched material from the 
interstellar medium (ISM; Alexander et al. 1998; Messenger 2000). 
For deuterium, the observed ratios are consistent with models and 
observations of the ISM, where low-temperature ion-molecule re- 
actions lead to enhanced D/H ratios in both gas- and solid-phase 
species (e.g. Millar, Bennett & Herbst 1989; Charnley, Tielens & 
Rodgers 1997.) However, there is little observational data on nitro- 
gen isotope ratios in the ISM, and models of the 15 N fractionation 
in typical dense clouds predict modest enhancements of ~ 25 per 
cent (Terzieva & Herbst 2000). In comparison, the largest 15 N en- 
hancements detected in meteorites - in so-called 'hotspots' - have 
15 N/ 14 N enhancements of more than a factor of four relative to the 
Earth (i.e. 8 15 N > 3000 %<0; Busemann et al. 2006). Values of 



<5 15 N(X) = 1000 x 
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where 



( 15 N/ 14 N) ffi is the terrestrial isotopic ratio of 0.003678 (De Bievre et al. 
1996) 



5 15 N > 1000 %o have also been found in hotspots in IDPs and 
Stardust samples (Floss et al. 2004, 2006; McKeegan et al. 2006). 

In an earlier paper we demonstrated that significantly in- 
creased 15 N fractionation can occur when CO is depleted onto dust 
grains (Charnley & Rodgers 2002 [Paper I]). The key fractionation 
reactions are: 



15 N+ 14 N 2 H+ 
15 N+ 14 N 2 H+ 



14 N+ 15 N 14 NH 4 



l N- 



N NH 4 
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where the exothermicities are AEi = 27.7 K and AE2 = 36.1 K 
(Terzieva & Herbst 2000). These preferentially drive 1 N into 
molecular nitrogen, at the expense of atomic N° which becomes 
isotopically light. In a standard dark cloud model the degree of 
fractionation is limited by chemical reactions which cycle nitro- 
gen between atomic and molecular form (see Fig. QJ. However, if 
CO is frozen out as ice, OH is unavailable, this cycle is broken, 
and much larger lo N-enhancements are possible. In this scenario, 
degradation of the ( 15 N-enriched) N2 by He + is essentially a one- 
way process. The N° atoms released by this reaction remain in the 
gas, whereas the N + ions react rapidly with H2 in a sequence of 
reactions to form ammonium (Fig. [TJ. The NH|" ions then recom- 
bine with electrons, producing NH2 and NH3, which subsequently 
freeze out on the dust grains. If the solid-phase NH2 radicals are 
hydrogenated to NH3, as expected from models of grain surface 
chemistry (Brown & Charnley 1990), the end result is a large abun- 
dance of isotopically heavy ammonia ice. At late times in this 
model, when N°/N2 > 1, the efficiency of reactions (TJ and O 
are greatly increased, leading to much larger gas-phase 15 N/ 14 N 
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Figure 1. Principal gas-phase reactions involving nitrogen. The dotted lines 
indicate the fractionation reactions {TJ and 



ratios than in the standard dense cloud chemistry. Ultimately, we 
found that the ammonia ice is enhanced in lj N by a factor of 1.8. 

This model could account for the presence of high bulk <5 15 N 
values in IDPs (Messenger 2000), and was also able to account 
for the non-detection of N2 in comets. Some additional process- 
ing was indicated to incorporate the 15 N-enriched ammonia into 
the carbonaceous matter, and secondary energetic processing to 
add lj Nfi2 side-groups to polycyclic aromatic hydrocarbon (PAH) 
molecules was suggested. Recent meteoritic studies have confirmed 
that most of the 15 N enrichment in IDPs is carried by amine side- 
groups on aromatic moeties (Keller et al. 2004). However, the 
model was unable to reproduce the largest <5 15 N values seen in the 
hotspots. Moreover, the experimental determination by Geppert et 
al. (2004) that recombination of N2H 4 " preferentially breaks the 
N=N bond, producing N and NH, acts to suppress the maximum 
fractionation (Rodgers & Charnley 2004 [Paper II]). 

In this paper, we are motivated to revisit our previous models 
as a result of several recent discoveries. Firstly, the experimental 
results of Geppert et al. (2004) have been challenged by Molek et 
al. (2007), who demonstrated that recombination of N2H + leads 
predominantly to N2 + H, as indicated in previous work (Adams 
et al. 1991). The new results indicate that rupture of the N2 bond 
occurs rarely, if at all, with an upper limit for the NH + N° branch- 
ing ratio of 5 per cent. Secondly, observations of several pre-stellar 
cores have revealed that they have temperatures below 10 K in their 
central, densest regions, as expected from physical models (Evans 
et al. 2001; Zucconi, Walmsley & Galli 2001). For example, Crapsi 



et al. (2007) derived T = 5.5 K in L1544, and Pagani et al. (2007) 
found T — 7 K in L134N. Due to the small zero-point energy 
changes associated with 15 N fractionation, reactions (TJ and l[2j are 
extremely sensitive to temperature, and these lower temperatures 
may be expected to yield larger 15 N/ 14 N ratios. 

Thirdly, observations of N2H + in dark clouds imply N2 abun- 
dances of a few x 10~ 6 (Womack, Ziurys & Wyckoff 1992; Maret, 
Bergin & Lada 2006). This is significantly less than the galac- 
tic elemental nitrogen abundance and Maret et al. proposed that 
the 'missing' nitrogen was in atomic form. Although the chemical 
model of Maret et al. assumed the Geppert et al. (2004) branching 
ratios for N2H + recombination, models assuming the 'old' branch- 
ing ratio also predict N° > N2 at early times (e.g. Pineau des 
Forets, Roueff & Flower 1990). This is because N — * N2 con- 
version occurs via neutral-neutral reactions, and the time-scale for 
this process, tn 2 , is larger than the dynamical free-fall time, ts. 
In contrast, CO is formed fairly rapidly by ion-molecule chem- 
istry, on a time-scale tco, where tco < tg < £n 2 - Therefore, 
we would expect young clouds to have high CO abundances but 
relatively small N2 /N ratios. This is consistent with the view that 
molecular clouds and pre-stellar cores form rapidly, on a time-scale 
of order a few times the free-fall time (Ballesteros-Paredes et al. 
2007). Dynamical-chemical models of collapsing clouds show that 
N — > N2 conversion is never efficient over tg (Brown, Charnley & 
Millar 1988). As discussed above, a large N°/N2 ratio is necessary 
to produce substantial 15 N-fractionation, so if the nitrogen in dense 
clouds is in fact mainly atomic rather than molecular, this will have 
important consequences for the 15 N/ 14 N ratios. 

Finally, in our earlier work we only calculated the bulk isotope 
ratios in the ammonia ice. In reality, the grain mantles will have 
an 'onion-ring' -like structure consisting of sequentially accreted 
monolayers. Thus, temporal variations in the gas-phase 15 N/ 14 N 
ratios will be preserved as spatial gradients in the ices, with each 
layer recording the gas-phase ratio at the time it was accreted. In 
particular, the late-accreting, uppermost layers will be the most 
highly fractionated. As it is these layers which are likely to expe- 
rience the largest degree of subsequent processing, it is necessary 
to distinguish between the bulk isotope ratios in the ice as a whole, 
and those in specific monolayers. 



2 MODEL 

We utilize the same chemical model as in our previous work (Pa- 
pers I & II). The lj N-fractionation chemistry is based on Terzieva 
& Herbst (2000); a list of the reactions and the rate coefficients 
we have adopted appeared in Paper II. We consider a high density 
core (riH 2 = 5 x 10 6 cm -3 ) with temperatures of 5, 7, and 10 K, 
and a cosmic ray ionization rate of 5 x 10~ 17 s _1 . The elemental 
15 N/ 14 N ratio in each model is generally assumed to be 1/400, al- 
though the effects of using other values were also investigated. A 
key parameter is the branching ratio for the dissociative recombi- 
nation of N2H + : 

N 2 H + + e — > N 2 + H (1-/) (3) 
— » NH + N (/) (4) 

We considered values for / = 0.0, 0.02, and 0.05, based on the ex- 
perimental results of Molek et al. (2007). We use elemental abun- 
dances for C, O, and N of 140, 290, and 80 parts per million respec- 
tively, relative to hydrogen (Savage & Sembach 1996), and assume 
complete depletion of metals. All the carbon is initially in the form 
of CO, and the remaining oxygen is atomic. 
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As before, we assume that CO freezes out onto grains with a 
sticking coefficient of unity, whereas N and N2 remain in the gas. 
The reason for this selective depletion of CO versus N2 is contro- 
versial, as laboratory experiments reveal that the two species have 
similar binding energies (Oberg et al. 2005). Nevertheless, obser- 
vations of pre-stellar cores show ample evidence for the presence 
of N2-rich, CO-poor regions toward the centers of these objects 
(Bergin & Tafalla 2007). Due to their low polarizabilities, atoms 
should have the lowest binding energies for physisorption to icy 
dust grains. Reaction with an H atom forms a simple hydride that 
can stick more effectively due to hydrogen bonding with the sub- 
strate. Thus, we assume that O atoms are hydrogenated on the 
grains to form water ice. Conversely, recent laboratory experiments 
indicate a very low efficiency for the reaction of N° and H on ice 
surfaces (T Hiraoka, private communication), which supports our 
assumption that N atoms do not stick. 

We have adapted our model to track the 15 N/ 14 N ratios as suc- 
cessive monolayers (ML) of ammonia ice are accreted. A 'typical' 
interstellar grain of size 0.1 /im has ~ 10 6 surface sites, and an 
abundance relative to hydrogen of ~ 10~ 12 . Hence, one ML corre- 
sponds to a total solid phase abundance of 10 -6 , and our elemental 
N abundance corresponds to a maximum of 80 ML of ice, assuming 
that all the nitrogen freezes out as NH3 (or NH2). In fact, a large 
fraction of the nitrogen remains in the gas phase in the form of N°, 
and we typically find that we form ~ 30 ML of ammonia ice. 



3 RESULTS 

We began by investigating the importance of the branching ratio, 
/, in the recombination of N2H + . We found only very small differ- 
ences in the results of models with / equal to 5, 2, and per cent, 
and conclude that, as long as / < 0.05, the presence of channel 
© does not significantly affect the 15 N chemistry. We then varied 
the elemental nitrogen isotope ratio, using values for 14 N/ 15 N of 
800, 400, and 100. In every case, the 15 N enhancements are identi- 
cal, relative to the underlying ratio. Henceforth, we discuss models 
with / = 0.02 and 14 N/ 15 N = 400. 



3.1 Temperature dependence 

Figure [2] shows the abundances of the three main repositories of 
nitrogen - gaseous N° and N2, and solid NH3 - as a function of 
time, assuming that all of the nitrogen is initially present as N2. 
The overall nitrogen chemistry is the same at 10 K and 7 K: N2 
is degraded on a time-scale ~ Myr. At late times just over half of 
the total nitrogen is in the gas phase as atoms, and just under half 
has frozen out as ammonia ice. At 5 K, however, very little ice is 
formed, due to the endo-ergocity of the reaction 



N + + H 2 



NH + H 



(5) 



We have assumed an activation energy barrier of 85 K for this re- 
action, based on the value in the UMIST reaction rate database 
(Woodall et al. 2007). At 5 K, the reaction becomes so slow that 
radiative recombination becomes the dominant loss route for N + 
ions, and only small amounts of gas-phase NH2 and NH3 are pro- 
duced. In effect, the barrier for reaction {5]l sets a limit on the ni- 
trogen chemistry in that, as the temperature drops, it eventually be- 
comes too cold to produce ammonia. In extremely cold cores the 
chemistry simply transforms the initial gas-phase N2 into gas-phase 
N°. In the following, we therefore look at the fractionation in 7 K 
and 10 K gas. 
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Figure 2. Time evolution of the main nitrogen carriers when all the N is 
initially molecular. Solid, dashed and dotted lines correspond to T = 10, 
7, and 5 K respectively. 



Figure [5^ shows the fractionation in the most important 
molecules. The results for 10 K gas are similar to those shown in 
Paper I, whereas at 7 K significantly larger peak 15 N/ 14 N ratios 
are found in gas-phase NH3 and N 2 . At early times, the gas-phase 
ammonia formed at 7 K is actually depleted in 15 N. This is due to 
the reduced rate of reaction l(5j which leads to a much larger N + 
abundance than at 10 K. These N + ions undergo isotope-exchange 
reactions with N 2 (Freysinger et al. 1994): 



■ J N+ + 14 N 2 



+ AE 3 



(6) 



where AE 3 = 28.3 K (Terzieva & Herbst 2000). This reac- 
tion shuffles lj N back into N 2 , reducing the 15 N-fractionation in 
N + and thus ammonia. At late times, this effect is still suppress- 
ing the N + fractionation relative to N2, but the enormous 15 N- 
enhancements in N2 ensures that N + and NH3 are more fraction- 
ated than at 10 K. In terms of the bulk 15 N/ 14 N ratio in the ammo- 
nia ice, both the 10 K and 7 K models predict similar enhancements 
of ~ 1.8. However, as discussed earlier, the time evolution of the 
isotope ratios will be preserved in the layered structure of the ices. 
Even though the overall icy 15 N/ 14 N ratios are the same in both 
models, ammonia ices which accreted at 7 K will have a more het- 
erogeneous structure, with the innermost layers depleted in 1 N, 
but with the top layers enhanced by factors ~ 10. 



3.2 The initial N/N 2 ratio 

Significant fractionation in N 2 can only occur when N2 is not the 
dominant nitrogen carrier. In the models discussed so far, this oc- 
curs at late times when most of the N2 has been broken down by 
He + attack, on time-scales ~ Myr. However, if a substantial frac- 
tion of the nitrogen is initially atomic, enhanced 15 N/ 14 N ratios can 
be produced on much shorter time-scales. The N°/N 2 ratio at t — 
depends on the dynamical and chemical history of the cloud before 
the formation of the dense pre-stellar core, and can be taken as a 
free parameter in our model. To investigate the importance of this 
factor, we have considered three additional models, where the ini- 
tial fraction of N in molecular form is taken to be 0.6, 0.1, and 0. 
The latter case, where all of the nitrogen is initially atomic, is not 
particularly realistic, but was used to constrain the amount of N 2 
that can be synthesized from N° in the limited time available be- 
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Figure 3. Time evolution of the 15 N/ N enhancement ratios in key species at T = 10 K (solid lines) and T = 7 K (dashed lines). The initial fraction of 
nitrogen in molecular form is 100%, 63%, and 10% in panels (a), (b), and (c) respectively. 



fore CO etc. freeze out. We find that this produces only very small 
N2 abundances, « 2 per cent of the total nitrogen. This means that 
the available pool of N2 - roughly half of which will end up as 
ammonia ice - is essentially equal to the N2 abundance at t — 0. 

Figure [3J3 and [5J; shows the isotope chemistry when a large 
fraction of the original nitrogen is atomic. Qualitatively, the gas- 
phase N2 and NH3 fractionation is similar to that in fig. [3^, except 
that the peak 15 N/ 14 N ratios occur much earlier as the initial N2 
abundance is reduced. As before, we find that the upper layers of 
the ice should be enhanced in 15 N by an order of magnitude. In 
terms of the overall bulk isotope ratio in the ice, it is clear that much 
larger values result, with ( 15 NHa/ 1 NH3)i co = 7 when only ten 
per cent of the initial nitrogen is molecular. This is due to the fact 
that less ammonia ice is formed overall, so the highly-fractionated 
upper layers account for a greater proportion of the total ice. We 
find roughly the same peak ratios regardless of the initial N°/N2 
ratio, but in models with more N2 the bulk fractionation is diluted 
by the large number of monolayers which accrete at early times 
with essentially normal (or even reduced in the 7 K model) 15 N/ 14 N 
ratios. 

3.3 Isotope ratios in individual monolayers 

We have calculated the 15 N/ 14 N in successive ML as they accrete 
from the gas. For comparison with laboratory measurements, we 
have converted the isotopic ratios into 5-values relative to the ter- 
restrial 15 N/ 14 N ratio. In order to do this, we need to know the orig- 
inal 15 N/ 14 N ratio in the material from which the protosolar neb- 
ula (PSN) was formed. We assume a value of ( 15 N/ 14 N) PSN = 
0.0025, which is implied by three independent measurements: 
Jupiter's atmosphere, high-temperature nebular condensates, and 
the solar wind (Fouchet et al. 2004; Meibom et al. 2007; Kallen- 
bach, Bamert & Hilchenbach 2007.) Figure [4] shows the 5-values 
in different ML for several models. Clearly, the upper layers of the 
ice are the most highly fractionated, with peak values of <5 15 N > 
3000 %o in every case. Also apparent is the approximate propor- 
tionality between the initial N2 abundance and the total number 
of ML that accrete. When only ten per cent of the nitrogen is ini- 
tially molecular, only four ML of NH3 ice are formed. However, 
in this case every ML is very highly fractionated, with <5 15 N every- 
where greater than 2000 %o. As discussed in £13-11 the peak values 



are larger for T = 7 K, but the 7 K model also produces a large 
number of ML at early times that are slightly depleted in 15 N. 



4 CONCLUSIONS 

We have revisited our earlier work on 15 N fractionation based on 
recent experimental and observation results. Using the branching 
ratios for N2H + dissociative recombination derived by Molek et 
al. (2007), we effectively recover the results of our earlier work in 
dense, CO-depleted gas at 10 K (Paper I.) We have investigated 
the effects of varying the temperature, and show that at lower tem- 
peratures, larger 15 N/ 14 N ratios are produced in gas-phase in N2. 
However, the barrier for the reaction of N + ions with H2 sets a 
lower limit on the temperature at which ammonia can be produced 
efficiently. Assuming the 'standard' rate coefficient for this reac- 
tion we find that very little ammonia ice is generated for T < 7 K. 
We have also looked at the effects of a substantial N°/N2 ratio at 
t = 0. We find that, because roughly half of the initial N2 ends up in 
the form of NH3 and NH2, reduced molecular nitrogen abundances 
yield less ammonia ice in total. Smaller N°/N2 ratios do not signif- 
icantly affect the peak gas-phase fractionation ratios, but because 
the highly-fractionated ammonia formed at late times represents a 
greater proportion of the total ice, we find that the bulk ice 15 N/ 14 N 
ratio can be greatly increased. 

Following the 15 N/ 14 N ratios in individual monolayers as they 
accrete sequentially from the gas, we have shown that gas-phase 
temporal variations in isotopic ratios are preserved as spatial gra- 
dients in the layered ammonia ice. The uppermost layers which 
accrete at late times have the largest 15 N-enhancements, up to an 
order of magnitude with respect to the elemental 15 N/ 14 N ratio. 
Converting to 5-values to compare with laboratory determinations 
of the isotope ratios in primitive solar system materials, we derive 
peak values of S 5 N > 3000 %o in gas at 10 K, and values as 
large as <5 15 N ~ 10000 %o in models with lower temperatures 
and smaller N°/N2 ratios. This is more than sufficient to account 
for the largest measured ratios, and demonstrates that interstellar 
gas-phase chemistry is likely the ultimate source of cometary and 
meteoritic 15 N anomalies. 

This work was supported by NASA's Origins of Solar Systems Pro- 
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Figure 4. 5 15 N values in different monolayers, for T = 10 K and T = 7 K. Labels refer to the fraction of nitrogen initially assumed to be molecular. 



gram through NASA Ames cooperative agreement NNX07AO86A 
with the SETI Institute, and by the NASA Goddard Center for As- 
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